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ABSTRACT 

A  navigation  system  concept  is  described  that  utilizes  the  moon  as  a  reflector 
in  a  bistatic  radar  system.  By  measuring  the  range  and  range  rate  of  the  receiver 
relative  to  the  moon,  the  location  of  the  receiver  on  the  earth's  surface  can  be  de¬ 
termined  in  both  longitude  and  latitude.  Lunar  radar  observations  have  shown  that 
lunar  range  and  range  rate  measurement  accuracies  equivalent  to  *30  m  can  be 
achieved.  By  placing  a  transponder  on  the  moon,  the  basic  measurement  accuracy 
could  be  improved  by  a  factor  of  5.  However,  for  operational  measurements  where 
rapid  readout  is  required,  these  accuracies  would  probably  be  degraded  by  a  factor 
of  5  to  10.  Placing  three  transmitters  at  appropriate  locations  on  the  earth's  sur¬ 
face  will  provide  worldwide  coverage.  A  transmitter  with  an  average  power  of  2 
MW  and  a  transmitting  aperture  of  170  m  would  supply  a  S/N  ratio  sufficient  for 
reliable  position  determination  with  a  dipole  antenna  receiver.  One  possible  radar 
receiver  configuration  incorporating  both  a  search  and  track  mode  is  given. 

A  mathematical  analysis  of  the  coverage  and  the  effective  position  accuracy 
indicates  that  (a)  worldwide  coverage  is  available,  (b)  time  coverage  is  restricted 
to  50  percent  on  the  average,  but  the  time  distribution  of  the  coverage  varies  over 
a  monthly  period,  (c)  at  low  latitudes,  the  effective  location  accuracy  is  a  function 
of  the  moon's  declination,  and  (d)  the  optimum  accuracy  is  obtainedat  high  latitudes. 

While  restricted  to  some  extent  in  coverage  as  described  above,  and  more 
sensitive  to  systematic  errors  in  the  ship's  velocity,  the  lunar  radar  navigation 
system,  when  compared  with  other  worldwide  radio  navigation  systems  such  as 
Omega  and  Transit,  can  achieve  higher  accuracies,  is  less  vulnerable  to  jamming, 
and  can  also  provide  an  independent,  one-way  communication  channel.  With  im¬ 
proved  technological  developments  and  active  reflectors  on  the  moon,  position  ac¬ 
curacies  of  the  order  of  several  meters  appear  possible. 


PROBLEM  STATUS 

This  is  a  final  report  on  one  phase  of  the  problem,  work  is  continuing  on  the 
project. 


AUTHORIZATION 
NRL  Problem  A01-35 
Project  33404,  Task  4880 


Manuscript  submitted  October  4,  1968. 
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A  LUNAR  RADAR  NAVIGATION  CONCEPT 
|  Unclassified  Title] 


INTRODUCTION 

Many  different  methods  have  been  employed  to  determine  the  position  of  ships  on  the 
ocean  or  of  other  unknown  locations  on  the  earth's  surface.  To  achieve  the  requirements 
of  increasing  accuracy,  and  all-weather,  global,  and  24-hour  coverage,  several  radio 
navigation  methods  have  been  developed.  Among  these,  the  Transit  satellites  and  the 
Omega  system  have  been  the  most  promising  under  development.  However,  practical 
considerations  have  necessitated  some  compromise  between  the  various  requirements. 

In  particular,  propagation  uncertainties  in  Omega  and  orbit  uncertainties  of  the  Transit 
satellites  limit  the  potential  accuracy  that  can  be  achieved. 

The  lunar  navigation  system  to  be  described  in  this  report  can  overcome  some  of 
the  limitations  of  the  other  systems  and  should  improve  the  potential  accuracy  by  an 
order  of  magnitude.  This  is  achieved  primarily  by  measuring  the  position  of  the  moon 
simultaneously  with  the  position  of  the  unknown  location.  To  obtain  optimum  accuracies 
(±5  m),  a  lunar  transponder  may  be  required.  In  addition,  the  system  is  less  vulnerable 
to  jamming  and  provides  one-way  communication  between  the  master  station  and  the  un¬ 
known  location.  The  basic  disadvantage  of  the  system  is  that  the  moon  is  above  the  hori¬ 
zon  for  only  about  50  percent  of  the  time. 


BASIC  CONCEPT  OF  THE  LUNAR  NAVIGATION  SYSTEM 

The  lunar  navigation  system  utilizes  the  moon  as  the  reflector  in  a  bistatic  radar 
system  where  the  transmitter  is  at  a  known  location  and  the  receiver  determines  its  po¬ 
sition  on  the  earth's  surface  from  measurements  of  ltd  distance  and  velocity  relative  to 
the  moon.  The  transmitting  station  measures  simultaneously  its  own  distance  and  veloc¬ 
ity  relative  to  the  moon  and  corrects  the  transmitted  signal  for  variations  of  the  lunar 
topography  and  the  lunar  ephemeris  errors.  The  receiving  station  derives  its  position 
by  applying  the  known  lunar  ephemeris  to  the  range  and  velocity  measurements.  The  two 
measurements  provide  the  necessary  corrections  to  both  longitude  and  latitude  of  the  un¬ 
known  location,  To  provide  worldwide  coverage  with  some  overlap,  three  ground-based 
transmitting  stations  are  required,  separated  by  120°  in  longitude  and  placed  at  moderate 
latitudes.  The  system  is  worldwide  except  for  the  limitations  imposed  by  the  time- 
varying  geometrical  configuration  of  the  earth-moon  system,  which  is  discussed  in  the 
next  two  sections. 


BASIC  LIMITATIONS  OF  COVERAGE 

There  is  a  basic  limitation  of  coverage,  due  to  the  availability  of  the  moon  at  a  given 
time  for  a  given  location.  The  moon  is  available  at  any  given  location  only  when  it  is 
above  the  horizon.  The  altitude  of  the  moon  (h)  at  a  given  location  is  given  by 

sin  h  -  sin  u>  sin  8  I  cos  q>  cos  8  cos  y  ,  (1) 


where 
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cp  is  the  geocentric  latitude  of  the  receiver, 

S  is  the  declination  of  the  moon,  and 
>-  is  the  local  hour  angle. 

The  limits  of  visibility  of  the  moon  are  given  by  Eq.  (i;  for  h  o.  This  yields  the 
usable  local  hour  angle  for  a  given  latitude  and  lunar  declination  as  expressed  in  Eq.  (2): 

cos  y  -  -ton  q>  tan  S  .  (2) 

In  Fig.  1,  the  usable  hour  angle  per  day  is  plotted  as  a  function  of  latitude  and  different 
declinations  of  the  moon.  For  a  declination  of  0°,  the  limiting  hour  angle  is  ±90°  or 
12  hours  for  all  latitudes.  For  positive  declination,  the  available  observation  time  per 
day  increases  with  latitude  until  24-hour  coverage  iB  obtained  at  high  latitudes.  At  large 
negative  latitudes  for  the  same  day  the  available  observation  time  is  reduced.  The 
opposite  behavior  as  a  function  of  latitude  occurs  for  negative  declination.  As  the  moon 
changes  its  declination  from  -28°  to  +28°  in  two  weeks,  the  moon  at  latitudes  close  to  the 
poles  will  be  observable  for  at  least  two  out  of  every  four  weeks. 


-90*  -70*  -SO*  -<I9*  -30*  -19  *  0  19  *  30*' 99*  60*  79*  90* 

I.ATITUDE  1^1 


Fig.  1  -  Visibility  of  the  moon  vs  latitude 
and  the  moon's  declination 


ERROR  ANALYSIS 

The  locus  of  all  points  on  the  earth's  surface  which  are  equidistant  from  the  moon  at 
a  given  time  is  a  circle  on  the  earth's  globe.  This  circle  is  located  in  a  plane  normal  to 
the  sublunar  distance  with  the  center  of  the  circle  positioned  below  the  sublunar  point. 
(See  Fig.  2.)  As  the  earth  rotates  around  Its  axis,  each  point  on  the  earth's  surface  will 
either  approach  or  recede  from  the  moon,  except  for  the  stationary  points  on  the  great 
circle  which  is  located  in  a  plane  passing  through  the  sublunar  point  and  the  earth's  axis 
of  rotation.  The  loci  of  equivelocity  points  are  small  circles  on  the  globe,  parallel  to  the 
plane  containing  the  great  circle  of  zero  velocity.  The  range  and  range  rate  measured 
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and  equivelocity  points 


TWO- LOCUTION  AMBIQUITY 
FOR  A  3IVEN  FREOUENCY 
SHIFT  AND  TIME  DELAY 


by  the  observer  will  establish  an  equidistance  and  an  equivelocity  circle.  The  intersec¬ 
tion  of  these  two  circles  determines  the  location  of  the  observer. 

In  general  the  intersection  of  the  equidistant  and  equivelocity  circles  corresponding 
to  one  measurement  will  provide  two  possible  locations  on  the  earth's  surface  (Fig.  2). 
However,  this  ambiguity  can  be  removed,  since  the  approximate  location  of  the  measured 
position  is  generally  known.  The  angle  of  intersection  of  the  two  circles  determines  the 
effective  accuracy  with  which  the  two  required  quantities  (longitude  and  latitude)  can  be 
determined.  In  general,  the  accuracy  of  the  longitude  is  less  dependent  than  that  of  the 
latitude  on  the  angle  of  intersection  of  the  two  circles.  The  latitude  error  Increases  as 
the  two  circles  approach  each  other  tangentially.  Since  the  sublunar  point  on  the  earth's 
surface  changes  with  both  a  daily  and  a  monthly  period,  the  relative  position  of  the  two 
circles  will  change  with  time,  and  thus  the  effective  error  of  the  longitude  and  the  lati¬ 
tude  will  be  a  function  of  both  time  and  position. 

In  Appendix  A,  the  effective  error  for  both  longitude  and  latitude  has  been  derived. 

In  the  computations  it  nas  been  assumed  that  the  measurement  errors  of  time  delay  and 
frequency  are  equal  for  the  most  efficient  operation,  and  that  the  corrections  for  errors 
in  the  lunar  ephemeris  and  for  variations  of  the  lunar  topography  have  been  made  by  the 
transmitting  station.  The  effective  errors  are  given  by  Eqs.  (3)  and  (4): 


"A*  = 


sin  <p  -  cos  ip  tun 


(3) 


a  \ "  — - — ■ — . . a ...  . 

''V  sin  tp  -  cos  ton  S  cos  y  m 

where  a^,  n\y,  and  am  are  the  standard  deviations  in  linear  dimensions  (meters)  of 
longitude,  latitude,  and  of  the  measurements,  respectively. 


These  expressions  give  the  effective  errorB  of  Ax  and  Ay  as  a  function  of  latitude  o, 
the  local  hour  angle  y,  and  the  moon's  declination  s.  For  some  limiting  cases  these 
errors  can  be  evaluated  in  a  simple  manner.  For  example,  for  n  =  o  and  --  la, 
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The  effective  error  of  the  longitude  is  constant  and  equals  the  measurement  accuracy  at 
all  locations.  Howeyer,  the  error  in  the  latitude  increases  as  the  latitude  approaches 
zero  (equator).  Similarly,  if  q>  o  and  fi  -t  o, 

1  m  __  1 _ m 

COS  )'  y  tan  ft  COS  y 

In  this  case,  the  errors  of  both  longitude  and  latitude  are  a  function  of  the  local  hour 
angle  and  are  minimum  at  transit. 

Data  for  a  complete  coverage  diagram  based  on  Eqs,  (3)  and  (4)  and  the  visibility  of 
the  moon  discussed  previously  are  der  ived  in  Appendix  A  and  given  in  Tables  A1  through 
A7,  In  general,  the  effective  accuracies  are  best  at  high  latitudes. 


MEASUREMENT  ERRORS 

The  effects  of  the  ionosphere  and  the  atmosphere  on  the  time  delay  and  frequency 
are  negligible  at  microwave  frequencies  for  measurement  errors  of  the  order  of  100  m, 
For  more  accurate  measurements,  corrections  for  typical  atmospheric  conditions  can 
reduce  errors  by  an  order  of  magnitude.  This  may  become  necessary  if  the  full  meas¬ 
urement  accuracy  of  the  lunar  transponder  is  to  be  utilized. 

The  main  modification  of  the  radar  signal  Is  caused  by  the  complex  reflection  char¬ 
acteristics  of  the  moon.  The  varying  topography  will  spread  the  leading  edge  by  extend¬ 
ing  the  buildup  of  the  pulse  anywhere  between  5  and  20  microseconds  (Fig.  3  shows  a 
typical  radar  echo).  In  addition,  the  frequency  is  spread  by  the  apparent  rotation  of  the 
moon  by  about  one  part  in  10 8.  However,  by  measuring  the  frequency  spread  in  the  first 
intervals  of  range  (about  10  to  100  psec),  the  spread  is  reduced  by  an  order  of  magnitude. 
The  measurement  consists  of  finding  and  establishing  the  centroid  of  both  the  leading 
edge  and  the  frequency  spectrum  with  an  accuracy  corresponding  to  1  fisec  in  time  and  a 
frequency  shift  of  one  part  in  1010.  This  is  probably  best  accomplished  by  utilizing  a 
tracking  gate  in  range  and  a  direct  frequency  count  to  establish  the  mean  time  delay  and 
frequency  of  the  received  signal. 


Fig.  3  -  Lunar  r$dar  echo  of 
May  15,  1968  (Observation  187, 
No.  75),  integration  time  0.1  sec 


CONFIDENTIAL 


CONFIDENTIAL 


NR  I,  UK  PORT  C814 


5 


Measurements  of  the  lunar  distance  with  the  NRL  lunar  radar  system,  which  em¬ 
ploys  a  basic  1.2-psec  range  resolution,  have  shown  that  the  internal  consistency  can  be 
maintained  to  within  tO.l  /isec  or  t30  m  in  flat  areas.  The  range  variations  for  a  lew 
adjacent  areas  are  given  in  Table  1.  However,  these  measurements  are  obtained  from 
the  fix  ft l  puini  uf  1  eileciiiMi  of  tiie  eciiu.  measuring  Lu  liie  mean  height  ot  a  lunar  area 
should  improve  the  consistency  of  the  range  measurements  and  reduce  the  frequency  of 
lunar  range  corrections  imposed  on  the  transmitter. 


Table  1 

Range  Variations  of  Several  Locations  on  the  Moon 


Time 

Interval 

L 

B 

V 

(m) 

General  Area 
on  Moon 

1/2  hour 

-6.95 

-7.11 

-1.72 

-1.72 

40 

Mosting 

1/2  hour 

-7.57 

-7.63 

-3.08 

-3.04 

-30 

Mosting 

1/2  hour 

6.85 

6.60 

5.73 

5.91 

-10 

Hyginus 

1  day 

1.43 

1.44 

6.78 

6.69 

-20 

Ukert 

1  month 

-3.79 

-3.78 

1.90 

1.71 

0 

Sinus  Medii 

Frequency  measurements  have  been  achieved  by  MIT  (1)  with  accuracies  of  one  part 
in  10 11  on  the  moon,  while  JPL  (2)  has  obtained  accuracies  of  the  order  of  one  part  in 
10 12  in  measurements  of  frequency  shift  with  transponders  from  the  Lunar  Orblter  se¬ 
ries.  The  frequency  spectrum  of  a  lunar  echo,  as  derived  from  an  NRL  cw  lunar  radar 
measurement,  is  shown  in  Fig.  4.  If  the  spectrum  from  only  the  Initial  return  is  ob¬ 
tained,  the  deviations  from  a  gausslan  spectrum  and  the  width  will  be  reduced,  thus  Im¬ 
proving  the  accuracy  with  which  the  centroid  of  the  spectrum  can  be  obtained. 


EFFECT  OF  SHIP'S  VELOCITY  ERROR  ON  POSITION  LOCATION 

Estimates  of  the  ship's  velocity  may  contain  either  random  errors  that  vary  with 
time  or  systematic  errors  that  remain  constant  for  several  hours.  Random  errors  with 
short  periods  (less  than  the  observational  Interval)  will  spread  the  frequency  spectrum 
and,  for  random  velocities  of  less  than  1  krn/hr,  will  not  materially  affect  the  accuracy 
of  the  position  determination.  The  effect  of  slowly  drifting  random  errors  of  the  ship's 
velocity  can  be  reduced  by  obtaining  several  measurements  in  order  to  average  the  effec¬ 
tive  position  error. 

The  systematic  errors  in  the  ship's  velocity  that  can  be  tolerated  without  degrading 
the  position  accuracy  are  proportional  to  the  effective  error,  as  discussed  in  the  error 
analysis  section. 
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FREQENCY  SPECTRUM  OP  A  LUNAR  RADAR  ECHO 


Fig.  4  -  Frequency  spectrum 
of  a  lunar  radar  echo 


The  effect  of  the  ship's  velocity  error  on  the  position  accuracy  has  been  derived  in 
Appendix  B  as  a  function  of  the  ratio  of  the  velocity  error  to  the  measurement  error. 
The  results  are  given  by 


8 


i 


"  \x  v 

"vT 


/il2(sinsj.  +  n^)1  1  "Av/u, 


and 


R2~-7~  siny(Sin3y  +  «JJ)*/J 

’Ay  « 


(b) 


(6) 


where 


°n 

aS3 

"Axy’  "Ay„ 
"Ax  •  "Ay 

m  m 

aAv/o)e'  °m 


the  local  hour  angle, 

cos  y  sin  o  -  cos  <p  tan  8, 

-  sin  y  sin  qs> 

the  errors  caused  by  the  velocity  error  on  the  longitude  and  latitude 
determination, 

the  errors  caused  by  the  measurement  errors  on  the  longitude  and  lati¬ 
tude  determination,  and 

the  velocity  and  measurement  errors  respectively. 
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Tables  B1  through  B7  represent  and  R  as  a  function  of  the  local  hour  angle  j  , 
the  geocentric  latitude  and  the  moon’s  declination  for  the  case  when  p,  Sv  )  •>„,  1. 

The  total  effective  errors  and  n  w  of  longitude  and  latitude  are  given  by 

r  ’r 

"Ur  v/l  • 

and 


"\v 


r 


/r  ■ » s 


For  velocity  errors  comparable  to  the  measurement  errors,  the  effective  position  error 
increases  on  the  average  by  about  20  percent.  For  the  assumed  measurement  error  of 
±300  m,  this  corresponds  to  a  velocity  error  of  ±3  cm/sec.  If  the  velocity  error  exceeds 
the  measurement  error  by  a  factor  of  four,  then  the  total  error  is  increased  by  100  per¬ 
cent  on  the  average.  For  larger  velocity  errors,  observations  over  several  hours  are 
required  to  determine  both  the  position  and  velocity  corrections,  if  the  full  measurement 
accuracy  is  desired. 


DESIGN  CONSIDERATIONS  FOR  LUNAR  RADAR 
NAVIGATION  SYSTEM 

In  considering  the  system  parameters  for  a  lunar  blstatic  radar  system,  it  has  been 
assumed  that  the  operational  measurement  accuracy  will  be  degraded  by  an  order  of 
magnitude  from  that  obtained  from  lunar  radar  astronomy  measurements  (1-3).  Fur¬ 
thermore  for  optimum  efficiency,  the  time  (range)  and  frequency  measurement  errors 
are  made  equal  as  discussed  In  the  section  on  measurement  errors.  This  corresponds 
to  a  1-psec  time  resolution  and  a  frequency  resolution  of  one  part  in  10 10,  which  are 
equivalent  to  a  position  error  of  ±300  m. 

Another  consideration  concerns  the  radar  frequency  to  be  employed  with  the  lunar 
radar  navigation  system.  Since  one  of  the  basic  measurements  is  the  doppler  shift  due  to 
the  relative  motion  of  the  ship's  location  with  respect  to  the  moon,  the  higher  the  radar 
frequency  the  greater  the  doppler  shift  and  in  turn  the  faster  a  readout  can  be  obtained 
for  a  given  accuracy.  On  the  other  hand,  the  complexity  of  a  transmitting  system  that 
supplies  sufficient  energy  to  a  dipole  receiving  system  increases  with  increasing  radar 
frequency.  These  considerations  will  require  a  compromise  between  the  complexity  of 
the  radar  transmitting  system  and  the  tolerable  observation  time.  With  larger  receiving 
antennas,  higher  frequencies,  or  more  complex  transmitters,  a  higher  accuracy  and  a 
more  rapid  readout  can  be  obtained.  An  indication  of  the  relation  between  these  factors 
and  frequency  is  given  by 


'A  t  r0  Ar 


where 

i  \i  =  the  measurement  accuracy  in  frequency, 
T0  =  the  observation  time, 

/ir  =  the  receiving  aperture, 
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30  m  with  an  average  power  of  60  kW.  This  will  also  reduce  the  required  readout  time 
by  a  factor  of  10. 

2.  The  time  at  the  transmitter  and  the  receiver  has  to  be  synchronized  within  the 
basic  measurement  accuracy.  Thus,  for  example,  for  an  accuracy  of  300  m,  the  two 
clocks  have  to  be  known  within  1  pee c.  The  use  of  atomic  clocks  is  clearly  indicated, 
and  periodic  resetting  of  the  receiver  clock  is  required.  For  a  clock  with  a  frequency 
stability  of  one  part  in  1013  (thallium),  the  receiver  clock  would  have  to  be  reset  every 
two  months. 

3.  Special  codes  may  be  used  to  reduce  the  time  of  search  in  the  uncertainty  interval 
of  the  range  measurements.  In  general,  a  search  should  be  completed  within  a  few  min¬ 
utes. 

4.  Special-purpose  computers  have  to  be  available  at  both  the  transmitting  station 
and  the  receiver  to  convert  the  known  information  of  the  lunar  ephemeris  and  the  position 
of  the  transmitter  to  range  and  range  rate  as  a  function  of  time.  The  expected  receiver 
position  is  introduced  to  determine  an  uncertainty  interval  from  which  the  corrected  po¬ 
sition  can  then  be  obtained. 

Operationally,  it  is  expected  that  the  transmitter  will  continuously  illuminate  the 
moon,  but  that  the  coding  interval  and  its  start  will  be  exactly  controlled  by  the  atomic 
clock  of  the  transmitter.  The  receiver  will  then  tune  In  on  the  received  signal  and  search 
for  its  position.  Once  the  position  is  found,  the  receiver  will  continue  to  track  the  posi¬ 
tion  and  indicate  the  motion  of  the  ship  or  vehicle  in  latitude  and  longitude.  The  effective 
error  of  the  Indicated  position  can  be  simultaneously  displayed  by  incorporating  in  the 
computer  the  information  presented  in  the  tables  of  Appendix  A.  Details  of  one  possible 
receiver  configuration  are  given  in  the  next  section. 


RECEIVING  SYSTEM 

The  ultimate  goal  in  designing  and  developing  the  receiving  system  is  a  package  that 
provides  automatic  latitude  and  longitude  information  with  little  or  no  human  adjustment. 
This  concept  is  possible  because  all  the  tradeoffs  for  receiver  simplicity  will  be  vested 
In  the  transmitter.  Since  simplicity  is  of  paramount  importance,  the  design  goal  for  the 
receiving  antenna  is  a  simple  dipole  or  whip  antenna.  This  eliminates  the  need  for  know¬ 
ing  where  to  point  the  antenna.  The  rest  of  the  receiving  system  will  utilize  Integrated 
circuits,  micromlnaturization,  and  modular  construction.  It  is  envisioned  that  the  final 
packaging  of  the  receiver  will  be  of  the  same  volume  and  size  as  today's  shipboard 
Loran  C  navigation  receivers. 

Figure  5  shows  a  proposed  block  diagram  of  the  receiving  system.  Since  biphase 
pulse  compression  is  assumed  in  the  transmitted  signal,  the  receiver  must  decode  the 
transmission  in  order  to  extract  the  necessary  information.  Also,  the  receiver  must 
pass  the  signal  through  fairly  narrow  filters,  which  requires  that  the  local  oscillator 
(LO)  chain  have  coarse  frequency  tracking  to  correct  for  the  earth's  rotation.  This  Is 
achieved  by  squaring  the  i-f  signal  to  remove  the  biphase  code  and  then,  with  a  frequency¬ 
following  filter  and  appropriate  counters,  developing  the  control  signals  for  adjusting  the 
LO  so  that  the  signal  will  pass  through  the  narrow  filters  used. 

Once  the  LO  frequency  is  correctly  set,  the  signal  is  fed  to  ten  correlators  and  fil¬ 
ters  (the  number  required  is  determined  by  operational  requirements)  to  determine  the 
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range  or  time  of  flight  of  the  signal.  The  code  generator  will  develop  the  proper  codes 
which  will  be  correlated  with  the  Incoming  signal  to  provide  coarse  range  information. 
The  range  tracker  must  sense  a  signal;  otherwise,  it  will  command  the  code  generator  to 
increment  the  code  until  a  signal  Is  found.  Once  the  signal  is  found,  the  range  tracker 
switches  into  the  fine  range  mode  to  measure  the  range  or  time  specified.  Also,  fine 
frequency  will  be  determined  after  the  fine  range  is  established,  and  this  Is  accomplished 
by  using  the  output  of  one  or  more  range  channels  and  feeding  its  signal  to  an  axis¬ 
crossing  counter  to  determine  precisely  the  frequency  of  the  received  signal.  The  range 
and  frequency  will  be  fed  into  a  computer,  which  will  display  latitude  and  longitude  auto¬ 
matically.  The  acquisition  time  is  determined  by  the  required  range  and  frequency  res¬ 
olutions  and  the  initial  uncertainty  of  the  position  under  observation.  Assuming  1-psec 
range  precision  and  100-km  position  uncertainty,  the  acquisition  time  is 


Maximum  (sec) 

Average  (sec) 

Lock  In  coarse  frequency 

40 

20 

Lock  In  coarse  range 

60 

30 

Lock  In  fine  range 

100 

50 

Measure  frequency 

100 

100 

Total  time  for  acquisition 

300 

200 

Once  the  signal  has  been  acquired,  updating  in  range  is  fairly  quick  (10  sec).  However, 
to  determine  the  frequency  will  still  require  100  sec.  The  times  used  are  for  a  particular 
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set  of  parameters.  Tradeoffs  in  performance  and  system  complexity  will  alter  the  ac¬ 
quisition  and  updating  times  required. 


DISCUSSION  AND  CONCLUSIONS 

To  evaluate  the  performance  of  the  lunar  radar  navigation  system,  it  may  help  to 
compare  it  with  the  Transit  system,  which  is  similar  In  fundamental  concept.  The  dif¬ 
ferences  of  the  lunar  navigation  system  with  respect  to  Transit  can  be  summarized  as 
follows; 

1.  A  natural  satellite  is  employed  as  the  reference  point  in  the  sky.  This  obviates 
the  necessity  of  periodically  Injecting  artificial  satellites  for  this  purpose,  but  restricts 
to  some  extent  the  coverage,  as  a  function  of  time  and  position,  which  would  be  available 
with  a  sufficient  number  of  artificial  satellites. 

2.  The  orbit  of  the  moon  is  well  known,  and  any  minor  corrections  to  its  orbit  can 
be  made  by  the  transmitter,  which  continuously  monitors  the  motion  of  the  moon.  This 
makes  It  possible  to  increase  the  effective  accuracy  of  the  range  and  range  rate  data  as 
compared  with  Transit,  where  the  orbit  of  the  satellite  has  to  be  updated  and  then  pre¬ 
dicted  for  a  period  of  several  hours. 

3.  The  relatively  slow  rotation  of  the  earth  as  compared  to  the  rapid  motion  of  the 
satellite  requires  a  higher  frequency  stability  at  the  receiver  for  a  comparable  accuracy 
in  position  measurements.  This  fact  also  makes  the  position  determination  more  sensi¬ 
tive  to  errors  In  the  ship's  velocity. 

4.  One-way  communication  is  available  with  some  minor  modification. 

5.  The  moon  is  not  vulnerable  to  hostile  Interference  and  is  much  more  difficult  to 

jam. 


6.  Finally,  due  to  the  large  distance  of  the  moon  from  the  earth,  the  system  may  be 
considered  In  the  future  as  a  navigational  aid  to  lunar  and  close  interplanetary  space¬ 
crafts. 

At  present,  the  achievable  accuracy  is  limited  technologically  by  the  accuracy  of  the 
available  atomic  clocks  and  the  need  to  reset  them  periodically.  With  the  future  devel¬ 
opment  of  atomic  clocks  and  the  placement  of  transponders  on  the  moon,  accuracies  of 
several  meters  should  be  possible.  On  the  other  hand,  for  stationary  locations  where  ob¬ 
servations  can  be  made  for  several  hours,  atomic  clocks  can  be  reset  automatically  and 
simultaneously  with  the  position  measurements. 

To  evaluate  more  closely  this  concept,  it  will  be  necessary  to  utilize  one  of  the  many 
radar  transmitters  now  available  for  testing  the  search  and  track  mode  of  a  receiver  In 
an  assumed  unknown  location  and  to  establish  more  definite  performance  criteria.  If  this 
concept  should  prove  Itself  practical,  it  may  then  be  feasible  to  consider  placing  a  trans¬ 
ponder  or  a  corner  reflector  on  the  moon  to  Increase  the  accuracy  as  well  as  reduce  the 
time  of  observation. 

Future  Improvements  of  the  lunar  radar  navigation  concept  may  include  the  utiliza¬ 
tion  of  additional  radio  reference  points  in  the  sky.  The  additional  radio  reference  points 
can  be  made  available  by  placing  isotropic  reflectors  at  the  llbration  points  of  the  moon. 
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The  applicable  libration  points  are  located  in  the  moon's  orbital  plane  and  consist  of  the 
two  stable  equilibrium  points  i-4  and  Ls  (4)  (equilateral  points)  at  ±60”  in  right  ascension 
relative  to  the  moon,  and  one  semistable  equilibrium  point  L ,  (collinear  nntnt)  opposite 
the  muon.  If  a  reflector  is  injected  at  any  one  of  these  points,  it  is  expected  to  stay  in¬ 
definitely  in  the  general  vicinity  of  the  particular  equilibrium  point.  The  small  periodic 
motions  of  the  reflectors  around  the  equilibrium  points  ran  he  compensated  for  by  the 
transmitter  corrections  as  before.  Utilization  of  these  reference  points  could  provide 
(a)  independent  determination  of  position  and  velocity  and  reduction  of  the  geometrical 
dilution  of  the  measurement  errors,  (b)  measurement  of  the  three  orthogonal  components 
in  both  position  and  velocity  for  limited  coverage,  and  (c)  an  increase  of  coverage  to  24 
hours  for  two  coordinates.  The  ability  to  determine  three  coordinates  in  both  position 
and  velocity  may  be  of  particular  use  when  the  concept  is  applied  as  an  independent  navi¬ 
gation  system  for  surface,  air,  or  spaceborne  vehicles,  but  the  extent  of  coverage  still 
has  to  be  investigated. 
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Fig.  Al  -  Geometry  of  the  earth-moon  system 


Ap  =  -  r  cos  9  cos  8  sin  y  A\ 

-  r  (cos  9  sin  S  -  sin  ip  cos  6  cos  y  )  A® 


(A3) 


and 


Ap  =  -  r  w,  cos  <p  cos  A'  cos  y  \k 

-  r  tu,  sin  <p  cos  5  sin  y  Aqj  .  (A4) 

Converting  the  angular  errors  to  linear  errors  by  setting  Ax  .  r  cos  <p  aa  and  Ay  =  r  a<p, 
the  errors  in  position  Ax  and  Ay  are  related  to  the  differentials  A p  and  A p  by  Eqs.  (A5) 
and  (A6): 


Ap 

cos  8 


sin  y  Ax  +  (cot  y  sin  9  -  cos  w  tan  8)  Ay 


(A5) 


Ap 

con  5 


-  coi  y  Ax  -  sin  y  sin  <j>  Ay 


(A6) 


It  is  seen  that  the  effective  measurement  error  is  a  weak  function  of  the  moon's  declina¬ 
tion  and  for  extreme  values  of  8  =  ±28°  will  Increase  the  basic  measurement  error  by 
about  10  percent.  To  obtain  the  corrections  of  Ax  and  Ay,  we  have  to  solve  Eqs.  (A5)  and 
(A8).  Let 
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and 


Therefore, 


t  cos3  q>  tan3  S 

-  2  sin'®  cos  ®  tan  £  cos  y)1/2. 


-  cos  y  cos  qa  tan  i 


°Ay  “  sin  9  -  cos  y  cos  q>  tan  £  ' 


These  expressions  will  give  .the.  effective  erroreof  sx  and- Ayas  a  function  Of  the  latitude 
e,  the  local  hour  angle  y ,  and  the  moon's  declination  a .  The  computer  printout  in  this 
appendix  gives  the  program  for  solving  Eqs,  (A13)  and  (A14). 


The  equations  are  computed  for  various  values  of  s ,  <p ,  and  y  in  Tables  A1  through 
A7  for  a„  s  i/cos  s.  The  tables  marked  "a"  indicate  the  effective  error  in  longitude, 
and  those  marked  "b"  give  the  corresponding  error  in  latitude.  The  first  vertical  column 
gives  the  latitude  of  the  observed  position  from  4-90°  to  -90°  in  10°  Intervals.  The  top 
row  gives  the  hour  angle  in  10°  intervals  (40  minutes  of  time)  starting  at  transit.  The 
tables  extend  symmetrically  for  negative  hour  angles  (not  shown).  Each  number,  when 
multiplied  by  the  given  basic  measurement  error,  indicates  the  effective  error  in  meters. 
The  tables  are  separated  by  about  two  days  in  time,  corresponding  to  a  change  in  decli¬ 
nation  by  9°.  The  effective  errors  increase  as  the  latitude  of  the  observed  position  ap¬ 
proaches  the  declination  of  the  moon. 
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Program  lor  Solution  of  Eqs.  (A13)  and  (A14) 
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TN9.4A 
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Table  A2a 

Error  in  Longitude  for  a;  Declination  Angle  of  18 
ERR9R  IN  LBNG1TU0E  FOR  DEC  -  IS 
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Error  in  Latitude  for  a  Declination  Angle  of  18 
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Table  A6b 

Error  in  Latitude  for  a  Declination  Angle  of  -18 
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Appendix  B 

EFFECT  OF  SYSTEMATIC  ERRORS  IN  SHIP’S  VELOCITY 
ON  POSITION  DETERMINATION 


When  the  location  of  a  moving  ship  is  to  be  determined,  it  is  necessary  to  include  in 
the  computation  the  doppler  shift  and  position  shift  caused  by  the  ship's  motion  relative 
to  the  moon.  The  error  in  the  ship's  velocity,  if  excessive,  will  be  directly  reflected  in 
an  error  of  position,  However,  by  taking  measurements  over  several  hours,  it  is  possi¬ 
ble  to  reduce  the  errors  in  both  position  and  velocity  to  the  basic  measurement  error. 

To  show  the  required  computation  as  well  as  the  functional  relation  between  position 
and  velocity,  the  observation  equation  will  be  derived.  The  range  and  range  rate  includ¬ 
ing  the  ship's  velocity  are  given  as  a  function  of  time  for  an  assumed  stationary  moon  by 


_ _  t 

p(t)  »  /d1  +  r1  -  2  Jr  cos  0(  t )  +  (  vsr  dt 


(Bl) 


Kt) 


d  (cos  0  ) 
~r  c?t  + 


(B2) 


where  v,r  is  the  radial  velocity  of  the  ship  relative  to  the  moon  and  the  other  symbols 
are  as  defined  in  Appendix  A.  For  a  short  observation  interval,  the  coordinates  of  the 
moon  and  the  position  of  the  ship  are  effectively  stationary,  and  it  is  further  assumed 
that  the  values  of  p  and  p  are  referred  to  the  midpoint  of  the  observational  interval  ro, 
or  (f-  t0)/2  =  r0/2.  Then 


pr  /t  «  7sJ  +  rJ  -  2  sr  cos  0  +  v,f  TQ/2 
d  (cos  0)1 


’V* 


dt 


+  v. 


Tn/t 


(B3) 

(B4) 


First,  v(r  will  be  derived  in  Greenwich  equatorial  coordinates.  Assume  a  coordi¬ 
nate  system  with  the  x  axis  pointing  to  the  Greenwich  meridian  and  the  *  axis  to  the 
north  pole  as  shown  in  Fig.  Bl.  Any  point  5  on  the  earth's  surface  can  be  expressed  in 
terms  of  the  earth’s  radius  r,  the  latitude  qp ,  and  the  longitude  A  (a  is  defined  as  the 
negative  of  the  conventional  astronomical  longitude). 

The  moon's  position  can  be  similarly  described  where  the  moon's  distance  is  s,  the 
latitude  qp  *  is  expressed  as  the  declination  s ,  and  the  longitude  a'  is  given  by  the  differ¬ 
ence  between  the  right  ascension  a  and  the  Greenwich  sidereal  time  o.  Thus  a'  =  a  -  e. 

The  radial  velocity  is  obtained  from 


V 


(B5) 
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Fig.  B1  -  Geometry  of  the  earth- moon  system 
Including  celestial  coordinates 


where  v,  la  the  shlp'B  velocity  vector,  and  the  negative  sign  indicates  that  the  radial  ve¬ 
locity  is  measured  relative  to  the  center  of  the  moon. 

To  determine  the  vector  v.  In  equatorial  coordinates,  first  define  v,  in  a  local  coor 
dlnate  system  at  the  point  of  observation,  with  the  y"  axis  pointing  east  and  the  z"  axis 
pointing  north.  The  x"  axis  Is  normal  to  the  y"  and  z"  axes,  and  forms  a  right-handed 
coordinate  system.  The  direction  of  v,  is  measured  by  the  angle  {  from  east  to  north 
as  shown  in  Fig.  B2. 


Fig.  B2  -  Components  of 
ship's  velooity 


The  components  of  »,  are  then  given  by 

vx„  =  0  , 

V  z  V,  cos  £  , 
v,„  a  v,  sin  i  . 

To  convert  to  the  Greenwich  equatorial  coordinate  system,  two  rotations  are  re¬ 
quired.  First  rotate  about  the  y"  axis  by  <p.  This  rotation  brings  the  x"  axis  parallel 
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to  the  equatorial  plane  and  the  z"  axis  parallel  to  the  earth's  axis  of  rotation.  A  second 
rotation  by  -a  about  the  z'  axis  brings  the  x‘  axis  to  the  Greenwich  meridian  and  makes 
ihp  y‘  avia  parallel  to  the  >  axis  in  Tig.  Dl.  Using  me  iransiurmation  equation  as  given 
by  Brouwer  and  Clemence,*  we  obtain  from  the  first  rotation  about  y“  axis  by  © 

vxt  -  vxu  cos  <p  -  v tH  sin  <p  =  -  vz„  sin  <p  , 

V  /  -  V  n 

y  y  ' 

vgl  =  vx«  sin  <p  +  v9„  cos  <p  =  vx„  cos  q>  , 
and  from  the  second  rotation  about  z'  axis  by  -a 

Vx  =  vxt  cos  A  -  Vyt  sin  A  , 

vy  =  vx»  sin  A  +  Vy>  cos  A  , 

V  =  V  l 

v  M  v  *  1  1 

giving  finally  the  velocity  vector  in  the  Greenwich-equatorial  coordinate  system  as 

v,  -  -  i  (vy*  sin  A  +  v,»  sin  ©  cos  A) 

+  1  (vyii  cos  A  -  v,»  sin  ©  sin  A  ) 

+  k  v,»  cos  ©  (B6) 

in  termu  of  the  east  and  north  velocities  vy»  and  v,„. 

The  lunar  distance  vector  is 

•  =  s  [i  cos  (a  -  9)  cos  S  +  J  sin  (a  -  9)  cos  5  +  k  sin  fi]  .  (B7) 

The  velocity  of  the  ship  relative  to  the  moon  is  then  given  by 


=  cos  B  [sin  y  vy »  +  (sin  ©  cos  y  -  tan  B  cos  ©)  v,»]  ,  (B8) 

where  y  is  the  local  hour  angle  of  the  moon,  l.e.,  y  =  a  -  a  -  a  and  a  is  now  defined  by  as¬ 

tronomical  convention. 

The  range  and  range  rate  equations  can  then  be  expressed  in  Greenwich  equatorial 
coordinates  as 

p  -  Vs*  +  rJ  -  2  ar  cos  (3 
T 

+  cos  5  £«in  y  VyH  +  (sin  ©  cos  y  -  tan  5  cos  ©)  vf»]  (B9) 

*D.  Brouwer  and  G.  M.  Clemenoe,  "Methods  of  Celestial  Mechanics,"  New  York:Academio  Press, 

p.  40,  1961. 
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and 

h  »  cos  S  [r«,  cos  v  sin  y  +  sin  y  vy" 
+  (sin  tp  cos  y  ~  ton  S  cos  p)  ■ 

The  differential  corrections  are  obtained  from 


(T 


8  +  *  ^.AV  +  K  "V>"  ' 

a;  |aa  +  £a*  *  +  4v‘“  • 

The  eight  coefficients  are  given  to  first-order  accuracy  (assume  (V, 

0V,»  )/2r  «  1  and  (T0/2)  Ayy„,  (Tc/2)Av,»  <  "„)l 


£p_ 

ex 

dp 

<5vy„ 


-  -  t  cos  o  cos  5  sin  y 


s  r  cos  5  (sin  <p  cos  y  -  cos  ip  tan  S) 


dv,„ 


tL  Z  -  cos  8  (rcu,  cos  <p  cos  y  +  cos  y  vy„  -  sin  ip  sin  y  V  ) 


aJL  =  -  COS  5  [r«„  sin  9  sin  y  -  (cos  9  cos  y  +  tnn  6  sin  q>)  v, 
d<Q 


«■— —  a  »in  y  COB  5 
tfVy  M 

tL  z  cos  S  (sin  <e  cos  y  *  cos  tp  t«n  5 )  • 

3v, » 

After  converting  to  linear  dimensions, 

Ax  =  r  cos  9  AX  , 
Ay  =  r  Ai  , 

the  observation  equations  become 


(BIO) 


a  )/2r  and 


.] 
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cos  S 


sin  y  Ax  +  (sin  p  cos  y  -  cos  up  tan  S )  Ay 


Ap 


cos  8 


=  f“  cos  v  ( - - - +  +  sin  y  sin  p - 

[  r  \rw0  cos  p  j  '  ro}„  cos  p 


Ax 


T  v«»  1 

+  -  sin  p  sin  y  +  (cos  p  cos  y  +  tan  8  sin  p)  -pjj—l  Ay 


+  sin  y  -  +  (sin  p  cos  y  -  cos  p  tan  8) 


Av.  h 


For  the  moderate  velocities  of  ships  vy„/(roie  cos  p),  vt«/(r&»#  cos  p),  and  vl»/(rtu1 
and  by  rearranging  the  equations  for  solution  of  Ax  and  Ay,  we  obtain 


Ap 

cos  8 


-  -  sin  y  Ax  +  (cos  y  sin  p-  cos  p  tan  8)  Ay 


and 


Ap 


AVy" 


A  v,» 


-  sin  y - (sin  p  cos  y  -  cos  p  tan  8) 

w.  cos  8  '  w,  '  », 


Let 


-  cos  y  Ax  -  sin  y  sin  p  Ay 


n  -  ! 

"i  "  cos  8 

Ap 


s'  ± 


sin 


Avv  * 


Av,» 


y  -  (cos  y  sin  p  -  cos  p  tan  8  )  < 


and 


then 


an  =  -  sin  y  , 

ala  a  cos  y  sin  p  -  cot  p  tan  8  . 
aal  s  -  cos  y  , 
a aa  -  -  sin  y  sin  p  i 


Ax  = 


ata 
ma  aaa 


(hill) 

(B12) 

,)  <c  1. 

(B13) 

(B14) 
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TN5.4A 


Computer  Program  Used  to  Solve  Eqs,  (B19)  and  (B20) 


09/0*/6fl 

coverage 

D I  HENS  I  ON  SIGX<19.7,l9),SlGY(JV,7,l9>, I  PRINT (3) 

DIMENSION  Rl (19,7,19)  ,R2(19./,19>,R3(l9.7»19),R4(l»,?.l9) 

C9MW0N/l/SIGX,Rl,RZ 

C9MMSN/2/SIGY«fi3,R4 

type  integer  dcln.ha 

DfcC,PHI,GAMMA,  TaND, SINP.C0SP, cSSG, DELTA  .COSO, SIND 
NUMERATR 
SING 

*  8H<«**,t>X.  ,0,6HFi».J))  ) 

ERRBP  IN  LONGITUDE  F0R  DEC  **.14,/) 


TYPE  DSUfllE 
TYPE  DSUBlE 
TYPE  DSUBLE 
DATA  < 1PRINT 
F8RMAT  <*l 
F8RMAT  <*O*,5x,*HA(*0R-)  0 


900 

901 

•105  120 

902  FSRMAT  <•♦* 

903  P0RMAT 

904  FORMAT 

905  F9RMAT 
*•*1/2 

906  F0RMAT 
•BC  *  ,14,/) 

907  f 0RMAT  (*1 
P8RMAT  («1 


135 
15) 

(*1 
(  *0*  > 
(79H1 
F8R 
(60M1 


150 

error 

Rl 
DEC 
R2 


165 


lb  30 

mo*,/,  • 


45  60 

PHI  *•/> 


75 


90 


in  latitude 

A12*< (A12**2 

i {6, / ) 

SJN0*(A12**2 


F0R  DEC  • • » 1 4 , / ) 

S I NG* *2 ) / <  A22 * • 2  •  Al2**2)) 
S I NG**2 ) *  *1/2  FOR  D 


908 

110 


DEC  ** , J 4 , / ) 
DEC  »* , 1 4 , / ) 


20 


30 

50 


EFFECT  8F  V  SUB  Z  0N  DELTA  Y  FOR 
EFFECT  0F  V  SUB  Y  ON  DELTA  Y  FOR 
FORMAT  (*5X,*,I2,*(1X*} 

DTR  •  0,0174532925199432960 

00  100  Jal.7 

DCLN  ■  9*  <1* J )  ♦  27 

DEC  ■  DCLN*DTR 

COSD  •  DCQS ( DEC ) 

SIND  •  OS  IN ( DEC ) 

TAND  •  SI ND/CSSD 
DO  50  I  ■  1,19 
LAT  ■  10*(1  -  1)  ♦  90 
PHI  •  LAT *DTH 
S I NP  *  DSIN(PHI) 

C9SP  ■  DCSS(PHI) 

DO  30  K  »  1,13 
HA  ■  15* («  •  1) 

OAMMA  i  HA*0TR 
C®SQ  ■  DCOS(GAMMA) 

S  J  NO  «  DS I N( GAMMA ) 

A12  ■  ABSF (C®SG*  s j  NP  -  C«SP*TAN0> 

A22  ■  •SINP*SING 

IF ( S I NP*S I ND  ♦  C6SP*ceSD*C8SQ  , 0 T ,  0)  Qfi  T(  20 
1 1 GX  1 1 ,  J,K  )  ■  S 1  GY ( I ,  J, K )  ■  R 1 ( 1 ,  J,  M )  *  R2<I,J,K>  ■  r3(I,j,K)  > 
»  R4 ( 1 1 j ,  K )  *  0, 

GO  TA  30 

DELTA  *  SI NP  «  CSSP*C0SQ*TAND 
IFtDELTA  , GO.  0)  DELTA  *  1.0D-30 

Rl ( I , J,K  >  ■  A12»SQRTF(A12**2  ♦  S1NQ»*2)/S3RTF(A22»*2  ♦  A12**2) 
R2 ( t , J,K )  «SING«SGRTF<Al2**2  ♦  SIN0**Z> 

1 1 GX  (  1  <  J  ,K  )  *  SGRTFd.  ♦  R1(I.U,K)**2) 

S I  G  Y  (  I ,  U  >  K  )  *  SQRTF  ( 1 ,  ♦  R2U,J.K)»*2) 

CONTINUE 

CONTINUE 

PRINT  905,  DCLN 
PRINT  901 
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TN5.4A 

00135  ]  «  1,19 

LAT  ■  10* (1  -  I )  ♦  90 
PRINT  902,  LAT 
M  ■  0 

00  133  L*1.13 

K  «  W  •  L 

If(Rl(I,J,K)  ,NE.  0, )  Q0  T0  134 

133  K  ■  K  •  1 

134  M  ■  K 

BNCSOE  (8,110. I  PR  1  NT  C2  )  )  M 

PRINT  IPRINT,  (fil  < I . J.K) ,K»1,M> 

135  PRINT  904 

PRINT  906.  DCLN 

PRINT  901 

00235  I  ■  1.19 

LAT  ■  10»<1  -  1 )  ♦  90 

PRINT  902,  LAI 

M  ■  0 

00  233  L«l, 13 
K  ■  14  .  L 

ir<R2<  I  , J,K>  , Nfc ,  o.)  Q8  T0  234 

233  K  a  K  -  1 

234  M  ■  K 

BNC0DE  (8,110. I  PR  1  NT <2> >  H 

PRINT  IPRINT,  <R2  ( |  ,  J,K) , Kal ,H ) 

235  PRINT  904 

PRINT  900,  OCLN 

PRINT  901 

00  35  I  ■  1,19 

LAT  a  10«(1  -  I)  +  90 

PRINT  902,  LAT 

M  a  Q 

DO  33  Kal, 13 

IFlSIQXd  ,J,K>  .60.  0.)  08  TO  34 

33  H  a  M  *  1 

34  ir(M  .EO,  0)  OH  TO  39 
BNCODE  (0,110.  IPRINT<2>)  14 

PRINT  IPRINT,  (SlQXd ,J,K>,Kal,M) 

35  PRINT  904 
PRINT  903,  OCLN 
PRINT  901 

00  40  |a  1,19 
LAT  a  10.(1  .  I J  ♦  90 
PRINT  902,  LAI 
M  a  0 

DO  38  Kal, 13 

irisiovii  ,j,kj  .es.  o ,  >  oe  to  39 

38  M  a  M  ♦  1 

39  IMH  .EO.  0)  00  T 0  40 

ENCODE  ( 8 , 110 . 1  HR  I  NT  C  2 )  )  M 

PRINT  IPRINT.  ( S  I G  V  < I >J,K),K»1,M) 

<0  PRINT  904 
100  CONTINUE 
BND 
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Table  Bla 

Position  Error  in  the  x  Coordinate  Caused  by  the  Velocity  Error  for  s  =  27 
Ai2*{( A12»*2  ♦  SIN6**2)/U22»»2  ♦  Al2**2) >*»l/2  FBR  DEC  = 
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Position  Error  in  the  y  Coordinate  Caused  by  the  Velocity  Error  for 
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Table  B2a 

Position  Error  in  the  x  Coordinate  Caused  by  the  Velocity  Error  for  s  =  18° 

R1  •  A12*((A12**2  *  SING»«2)/(A22**2  *  A12**2> )«»J/2  F8R  Dec  *  18 

HA(*#R-i  0  15  30  45  60  /5  90  105  120  135  150  165  180 
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Table  B3b 

Position  Error  in  the  y  Coordinate  Caused  by  the  Velocity  Error  for  8  =  9° 

R2  *  SINR*(A12**2  ♦  SING**2>«»l/2  F$R  DEC  e  9 
HA«*iR-l  0  15  30  45  60  /5  90  105  120  135  ISO  165  180 
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Table  B4a 

PosiUon  Error  in  tbe  X  Coordinate  Caused  by  the  Velocity  Error  for  «=0° 

«1  •  A12*((A12**2  *  S1N6**2)/(A22«*2  ♦  A12**2)  >*U/2  F9R  DEC  .  fl 

«<•«-.  0  15  «  45  „  „  „  ,3,  1M  li( 
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Table  B4b 

Position  Error  in  the  y  Coordinate  Caused  by  the  Velocity  Error  for  s  = 
5lNG»( A12*»2  ♦  SIN6**2)**l/2  1 8  R  DEC  =  0 


CONFIDENTIAL 


NR I.  REPORT  6814 


49 


Table  B6a 

i  Error  In  the  *  Coordinate  Caused  by  the  Velocity  Error  for 
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Position  Error  In  the  y  Coordinate  Caused  by  the  Velocity  Error  far  8=  18 
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Position  Error  in  the  y  Coordinate  Caused  by  the  Velocity  Error  for 
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13  ABSTRACT 


A  navigation  system  concept  is  described  that  utilizes  the  moon  as  a  reflector  In  a  bistatic 
radar  system.  By  measuring  the  range  and  range  rate  of  the  receiver  relative  to  the  moon, 
the  looatton  of  the  receiver  on  the  earth's  surface  can  be  determined  in  both  longitude  and  lati¬ 
tude.  Lunar  radar  observations  have  shown  that  lunar  range  and  range  rate  measurement  ac¬ 
curacies  equivalent  to  ±30  m  can  be  aohieved.  By  placing  a  transponder  on  the  moon,  the  basic 
measurement  accuracy  could  be  improved  by  a  factor  of  6.  However,  for  operational  meas¬ 
urements  where  rapid  readout  is  required,  these  accuracies  would  probably  be  degraded  by  a 
factor  of  5  to  10.  r lacing  three  transmitters  at  appropriate  locations  on  the  earth's  surface 
will  provide  worldwide  coverage.  A  transmitter  wlla  an  average  power  of  2  MW  and  a  trans¬ 
mitting  aperture  of  170  m  would  supply  a  S/N  ratio  sufficient  for  reliable  position  determina¬ 
tion  with  a  dipole  antenna  reoeiver.  One  possible  radar  receiver  configuration  incorporating 
both  a  search  and  track  mode  Is  given. 

A  mathematical  analysis  of  the  coverage  and  the  effective  position  accuraoy  indicates  that 
(a)  worldwide  coverage  is  available,  (b)  time  coverage  is  restricted  to  60  percent  on  the  aver¬ 
age,  but  the  time  distribution  of  the  coverage  varies  over  a  monthly  period,  (c)  at  low  latitudes, 
the  effective  location  accuracy  is  a  function  of  the  moon's  declination,  and  (d)  the  optimum  ac¬ 
curacy  la  obtained  at  high  latitudes. 
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While  restricted  to  some  oxteut  in  coverage  as  described  above,  and  more  sensitive  to 
systematic  errors  In  the  ship's  velocity,  the  lunar  radar  navigation  system,  when  compared 
with  other  worldwide  radio  navigation  systems  such  as  Omega  and  Transit,  can  achieve  higher 
accuracies,  is  lese  vulnerable  to  Jamming,  and  can  also  provide  an  independent,  one-way 
communication  channel.  With  improved  technological  developments  and  active  reflectors  on 
the  moon,  position  accuracies  of  the  order  of  several  meters  appear  possible. 
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